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Abstract 

In cosmology minisuperspace models are described by nonlinear time-reparame- 
trization invariant systems with a finite number of degrees of freedom. Often these 
models are not explicitly integrable and cannot be quantized exactly. Having this 
in mind, we present a scheme for the (approximate) quantization of perturbed, non- 
integrable, time-reparametrization invariant systems that uses (approximate) gauge 
invariant quantities. We apply the scheme to a couple of simple quantum cosmological 
models. 
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1 Introduction 

Historically, the first attempt in quantizing gravity moved from the canonical formalism. The 
origin of quantum gravity dates back to the pioneeristic works of Bergmann and Goldberg [|l| 
and Dirac [0 who first investigated in the late 50 's the canonical analysis of classical general 
relativity and applied to the gravitational field the theory of quantization of constrained 
systems This line of research led to the Hamiltonian formulation of the Einstein theory 
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by Arnowitt, Deser, and Misner ^ and to the application of Dirac quantization to the ADM 
formahsm by Wheeler [^] and DeWitt (Wheeler-DeWitt equation). 

In spite of an effort lasting more than 40 years, the conceptual problems arising in the 
quantum canonical formulation of gravity are hale and hearty and a large amount of time 
is presently devoted to the discussion of "structural issues in quantum gravity" [[7|. Let us 
just mention the so-called "problem of time" p, |^: since the first-class Hamiltonian for 
general relativity is identically zero, due to general covariance, the evolution of the system 
is "hidden" in the constraints. 

Minisuperspace models are the natural arena to investigate these structural issues. By 
retaining only a finite number of degrees of freedom quantum field theory reduces to quantum 
mechanics and typical problems due to the field nature of the system (for instance, anomalies 
see Ref. |]10|) disappear. However, other conceptual problems, as the mentioned "problem 



of time", survive to the reduction to a finite number of degrees of freedom [0, |T2|. The 
standard techniques of quantization of constrained gauge systems cannot be naively applied, 
for instance, even to the simple case of a Friedmann-Robertson- Walker (FRW) universe 
minimally coupled to a massive scalar field: for nonintegrable systems the absence of a 
global time parameter forbids the removal of the gauge degree of freedom and a sensible 
quantum theory cannot be obtained. So some approximation must be introduced - at some 
stage - in the quantization procedure. A possible approach is given by the semiclassical 
approximation. In this framework a (approximate) Schrodinger equation for the system can 



be derived |]T3|. However, the definition of time in the Schrodinger equation - and the Schro- 
dinger equation itself - depend on the particular wave function considered. As a consequence, 
it is not obvious how to construct a Hilbert space. 

In this paper we present a method to obtain an approximate global time definition and 
approximate gauge invariant quantities for finite-dimensional, nonintegrable, perturbed sys- 
tems, i.e. systems whose nonintegrable sector can be written as a small perturbation. We 
develop a scheme that allows the explicit construction of these quantities to any order in the 
perturbation parameter A. Then the approximate system to the order 0(A") can be quan- 
tized using the standard operator approach. Finally, we illustrate the method discussing a 
couple of quantum cosmological models. 



2 Quantization of Integrable Finite-Dimensional Time- 
Reparametrization Invariant Systems 

As a warming up exercise let us discuss finite-dimensional systems that are explicitly inte- 
grable. Due to time-reparametrization invariance their action reads 



S= r dt [%p,-u{t)H{q,p)] , (2.1) 
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where /i = 0, l...A^ — 1 and u{t) is a nondynamical variable that imposes the constraint 
H = 0. It is straightforward to see that the action ( |2.1| ) remains invariant under the time 
redefinition t ^ i = fit) |T^. This property imphes that the system is invariant under the 
(gauge) transformation generated by H 

dH 

dH 

6p^ = -e— = e[p^,H]p, (2.2) 



^ de 
ou = — 
dt 



where e{ti) = 6(^2) = 0. Note that the gauge transformation (|2.2| ) has the same content of 
the equations of motion. Since the system is integrable Eqs. (|2.2|) can be exphcitly solved. 



The standard operator approaches to the quantization of the system ( |2.1|) are the Dirac 
method and the reduced method. In both approaches the redundant gauge degree of freedom 
is eliminated via the imposition of an extra condition F{q,p; t) = (gauge identity) [|15]: in 
the Dirac method the constraint H = is promoted to a quantum operator and then the 
gauge degree of freedom is eliminated by the gauge fixing; conversely, the reduced method 
leads first to a classical reduced gauge-fixed phase space where quantization can be carried 
out as usual (Schrodinger equation) and wave functions have the customary interpretation. 

Since the gauge fixing condition is given in the form of a local relation, it may happen 
that the gauge fixing condition does not globally hold in the entire phase space (Gribov 
obstruction, see for instance Ref. [Q)- So the quantization cannot be carried out. (Geomet- 
rically, the existence of a global time parameter implies that gauge conditions intersect the 
gauge orbits on the constraint surface once and only once.) However, when the system is 
integrable it is certainly possible to choose a global gauge fixing. Indeed, we can construct a 
canonical transformation {qfi,p^} — > {Q^, P^} where one of the new canonical coordinates, 
say Po, coincides with the constraint, its conjugate Qo transforms linearly under the gauge 
transformations, and the remaining canonical coordinates are gauge invariant (Shanmugad- 
hasan canonical variables) |T^. The conjugate variable to the Hamiltonian, Qq, defines a 
global time parameter and there is no Gribov obstruction. Let us see this in detail. 

Consider Eq. ( ^.1| ) and neglect the constraint H = 0. Since the system is integrable we 
can write and solve the equation for the Hamilton characteristic function 

Hiq^,^\=ao, ^ Vr = W^(g^,M. (2.3) 

where are N constants (of motion). The Hamilton characteristic function generates a 
canonical transformation to cyclic coordinates {Q^, P^} 

dW 

Qo = j^ = r + Po, (2.4) 
oao 
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where 



dW 

Q. = ^ = A, ^ = 1,2,... (2.5) 

= «M > (2-6) 

(2.7) 

r(t) = f u{t')dt' . (2.8) 

The quantities {Q^, P^} form a set of Shanmugadhasan variables of the system. Indeed, 
since are constants of motion, Qi and Pi are 2(A^ — 1) gauge-invariant quantities, Pq is 
the Hamiltonian and Qo transforms hnearly for the gauge transformation generated by H. 

and can be written as functions of q and p evaluating Eqs. ( |2.4| - |2l^ ) at r = and 
inverting the relations 

dW dW 
oaf, dq^ 

The canonical transformation to Shanmugadhasan variables finally reads 

Qf. = P^^{q,p), P^ = a^{q,p). (2.10) 

The quantities {Qi, P^} form a complete set of gauge-invariant quantities (observables) that 
are in a one-to-one correspondence with the initial conditions. Note that the definition of 
Shanmugadhasan variables given by Eq. ( |2.10 ) is not unique. For instance, the generator 



of the canonical transformation is defined up to an additive generic function of the new 
momenta f{a^) 

df 

P^^P'^ = P^, Q^^^Q'^ = Q^ + -L.. (2.11) 

Clearly, {Q'^, P'^} form a different set of Shanmugadhasan variables. 
In the Shanmugadhasan representation ( |2.10D the action reads 

S = [ dtiQoPo + Q^P^ - uPo) . (2.12) 



The quantity Qq can be used to fix the gauge because its transformation properties under 
the gauge transformation imply that time defined by this variable covers once and only once 
the symplectic manifold, i.e. time defined by Qo is a global time (see e.g. Ref. [|17|)- Using 
the Shanmugadhasan variables the quantization procedure becomes trivial and both Dirac 
and reduced approaches lead to the same Hilbert space. 

It is well-known that classical canonical transformations and canonical quantization gen- 
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erally do not commute. Graphically (operators are marked with the hat symbol) 



H(q,p) Classical 
Theory 



Quantization 
Algorithm 



H{q,p) Quantum 
Theory 



^ Canonical ^ 
Transformation 



??? 



H{Q,P) Classical 
Theory 



Quantization 
Algorithm 



H{Q,P) Quantum 
Theory 



(2.13) 



So the quantization of a classical system described by different - but classically equivalent - 
sets of canonical coordinates may lead to different quantum theories. The Shanmugadhasan 
variables are not immune from this disease: quantization in the {Q, P} representation may 
be inequivalent to quantization in the {q,p} representation. However, the Shanmugadhasan 
set of canonical variables seems to represent a preferred choice of canonical coordinates to 
be used in the quantization procedure. Indeed, using the Shanmugadhasan variables both 
Dirac and reduced approaches lead to the same Hilbert space. This property is not true in a 
generic representation and gives the Shanmugadhasan set of coordinates a preferred status. 



In the Dirac approach we start setting the commutation relations 



(2.14) 



The quantization of the system can be achieved imposing that the constraint becomes an 
operator on wave functions. This enforces the gauge invariance of the theory. Once one has 
imposed the quantum relation H"^ = 0, the gauge must be fixed. This can be done using 
the Faddeev-Popov procedure to define the inner product. We have 



(^2,^1 



d[Q,] mQ^)5iF)AFp^^iQ 



(2.15) 



where d[Qfj] is defined in the unconstrained phase space and represents the measure invariant 
under the symmetry transformations of the system (rigid and gauge transformations). App 
is the Faddeev-Popov determinant (see e.g. Ref. |T^ ). 



Using the measure d[Q^] = JJdQ^ (we neglect possible issues related to the self-adjoint- 
ness nature of the operators) and the operator representation Qfj., —id^ the 

constraint equation reads —ido'^i^Qi) = 0. So the physical wave functions do not depend 
on Qq and the zeroth degree of freedom is pure gauge. A suitable gauge fixing condition is 
F = Qq — k = 0, where is a parameter. With this choice the Faddeev-Popov determinant 
is AiT'p = 1 and the inner product ( ^.151 ) reads 



(2.16) 
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This completes the quantization of the system. The Hilbert space can also be obtained 
reducing first the system to the physical degrees of freedom by the gauge fixing and then 
applying the quantization algorithm. In this case we impose the gauge fixing condition 
F = Qo — t = 0. This determines the Lagrange multipher as m = 1 because of the gauge 
equations and of the definition of Qq. Now the gauge fixing can be implemented by the (time 
dependent) canonical transformation 

Q'o = Qo-t, P;, = Po, (2.17) 

imposing Qq = (gauge fixing) and Po=0 (constraint). The effective Hamiltonian on the 
gauge shell becomes H^s = Pq = and the Schrodinger equation reads 

^|^^(t,g,) =0. (2.18) 

Wave functions do not depend on t = Qo and we recover the same Hilbert space that we 
have previously found following the Dirac approach. 

We have seen that for integrable systems a maximal set of gauge invariant canonical 
variables can be constructed. Then the quantization is carried out using the new variables. 
Now let us discuss a concrete example. Consider the Hamiltonian constraint 

H = -lipl + ql) + l{pl + ql)=0. (2.19) 

(Equation ( |2.19| ) describes, for instance, a closed FRW universe coupled to a homogeneus 
conformal scalar field. In this case go is the scale factor of the metric a and x = O'l/cf is the 
conformal scalar field UTTl.) 



The Hamilton characteristic function can be chosen (we consider without loss of generality 
Po and pi positive defined) 

W = J dqo,j2{ai-ao)-qi + J dqi^2ai~qj. (2.20) 

Using Eqs. ( |2.9| ) it is straightforward to obtain the Shanmugadhasan variables 

Qo = arccos ^i^ = , (2.21) 

Ql = — arccos ^° — arccos , '^^ , (2.22) 

VPo + 9o ypi + Ql 

Po = H{q,,p,), (2.23) 

Pi = lipi + q!)- (2.24) 

The system can be quantized using Qq as time parameter. In the original variables the time 
parameter reads 

t = arctan — . (2.25) 
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In this case the system is separable and time is defined only by one of the two degrees 
of freedom. So we can reduce the system to the Shanmugadhasan form only in the sector 
defined by {qo,po)- In this representation the "hybrid" canonical variables are (Qo, Pq, qi,Pi), 
where 

Q'q = arccos . ^° , (2.26) 
VPo + Qq 

Po = -liPl + ^o)- (2.27) 
In the hybrid representation the Schrodinger equation reads 

t) = ^ (^-|^ + ql^ tlj{qi; t) , (2.28) 

where we have used the usual representation for the operators qi and pi. The eigenfunctions 
of the Hamiltonian are 

V^„(gi) = ^ H4qi)e-'^y^ , (2.29) 

where Hn are the Hermite polynomial of order n. The eigenfunctions (|2.29|) form a orthonor- 
mal basis in the Hilbert space. 



3 Quantization of Perturbed Nonintegrable Systems 

Any integrable time-reparametrization invariant system with a finite number of degrees of 
freedom can be exactly quantized in the Shanmugadhasan (or hybrid) form. The canonical 
coordinate conjugate to the Hamiltonian can be used to define the time parameter for the 
system and a Schrodinger equation can be written. However, the quantization procedure 
illustrated in the previous section relies deeply on integrability and on the existence of a global 
time. If the system under consideration does not admit a global time exact quantization is 
not possible. So an approximate method of quantization for these systems is worth exploring. 

In this paper we are interested in systems that are described by a Hamiltonian constraint 
of the form 

H = Hintiqi^.Pf,) + -^^pert('?M'PM) ' (3-1) 

where A is a small dimensionless parameter, i^int identifies the integrable sector of the theory, 
and -ffpert represents the nonintegrable perturbation. 



We have mentioned in the introduction the semiclassical approximation ||13|. In this 
paper we follow a different approach. Our aim is to obtain approximate gauge invariant 
quantities for the system, i.e. approximate Shanmugadhasan variables up to a certain order 
n in the parameter A. The action is reduced to the form (|2.12|) plus a term of the order 
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0{X'^~^^) and the (approximate) system can be quantized along the hnes of the previous 
section. Let us see the procedure in detail. 

Since the first term in the Hamiltonian is completely integrable, let us suppose without 
loss of generality that the Hamiltonian is 

H = po + XH,{q^,p^). (3.2) 

(Is it always possible to reduce Eq. ( p.l| ) in the form ( |3.2| ) using the canonical transformation 
described in the previous section.) The reduction to the Shanmugadhasan form is obtained 
implementing the canonical transformation generated by 

W{q^, P,) = q,P, + XW^{q,, P,) + \^W2{q^.. P,) + ...., (3.3) 

where Wi{q^, P^), W2{q^, P^i),--- are determined by the request that the transformed Hamil- 
tonian has the form 

H = Po + 0(A"+i) . (3.4) 

Calculating and from Eq. (p.3|), expanding in powers of A, and equating terms of the 
same order we find 

dWi 

+ Hi{q^,p^)\p^=P^ = 0, 

From Eqs. ( p.5|) it is straightforward to obtain the expressions for Wn- For instance Wi reads 

Wi{q^,P^) = - J d%Hi{qQ,qi,p^)\p^=p^ + Z{qi,P^) , (3.6) 

where Z is an arbitrary function that we set for simplicity equal to zero. (Note that all Wn are 
determined up to additive arbitrary functions of the new momenta: Wn Wn + Zn{P^) since 
Eqs. ( |3.5| ) involve only derivatives of Wn w.r.t. q.) The first-order canonical transformation 
is generated by 

W{q,,P,)=q,P,-Xj dq',H,{q',,q,,p^)\,^=p^ + 0{X'). (3.7) 
The second term reads 

«M.,p.)^f.,i{..^.r.,stf)_, (3,s) 

and so on. Implementing the canonical transformation generated by Eq. (|3.3| ) to the order 
n one finds approximate gauge invariant variables to the order ©(A"). The Poisson brackets 
read 

[Q^,P.]p = 5^u + 0iX^^'). (3.9) 
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Finally, the canonical variable conjugate to Pq 



truncated to the order 0(A") is the approximate time parameter of the system. 

We do not expect in general that the series ( p.lO| ) are converging in the entire phase space 
because a nonintegrable system does not admit a global time (see e.g. Ref. [0). Nevertheless 
Qo defined in Eq. ( p.lQ| ) is a good time parameter in those regions of the phase space where 
the Poisson bracket of Qo with the Hamiltonian is positive defined, i.e. when the quantity 
0(A"+^) in Eq. ( p. 91) is less than 1. This happens usually when the perturbation ifpert is 
small compared to the integrable part. Conversely, when iYpert ^ 1 nonperturbative effects 
are present and the series (|3.1CI|) will in general not converge. Clearly, for A = the condition 



is satisfied and the zero-order term of Qo in Eq. ( p.lO| ) is a global time, as expected because 
now the system is integrable. 

At this point one can quantize the system using the (approximate) Shanmugadhasan 
variables along the lines of Sect. |^. The quantization procedure can be completed exactly 
and formally does not depend on the expansion scheme. Indeed, the expansion is chosen 
such that the new Hamiltonian is Pq and the wave functions are diagonalized by suitable 
operators defined by the gauge invariant variables Qi and Pi. The approximation scheme af- 
fects only the definition of the new canonical quantities as functions of the original variables 
of the system. In conclusion, the procedure simply allows to calculate the finite (approxi- 
mate) canonical transformation to a set of Shanmugadhasan variables. For nonintegrable, 
time-reparametrization invariant, perturbed systems time and gauge invariant quantities are 
approximate concepts defined only in suitable regions of the phase space. 

In the next section we shall apply the method illustrated above to a simple minisuperspace 
model. 



4 A Flat FRW Universe Coupled to a Scalar Field 



Let us consider a fiat FRW universe coupled to a scalar field ip. In Planck units (/pianck 



AnG/S) the action density for this model reads 

10V 



S 



dt 



1 aa^ 



2 



(4.1) 



where a > is the scale factor of the FRW universe, N is the lapse function, V{ip) is the 
potential of the scalar field, and dots represent derivatives w.r.t. the time t. Equation (|4.1| ) 
can be cast in the Hamiltonian form. Defining u = N/ we have 



S = J dt japa + ^p^-u ^ (pI - o^pI) + aV(^) | 



(4.2) 
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Clearly, the action ( [4.2|) has the form (|2.1| ); u is the nondynamical variable. When the 
potential is constant the system is separable and integrable; conversely, when the potential 
for the scalar field is not constant the system is not integrable. However, if the conditions 

aV(^)«aV^^pJ (4.3) 

are satisfied, the last term in Eq. ( [4. 2] ) can be considered as a small perturbation and we can 
apply the techniques developed in the previous section. (For instance, a massive minimally 
coupled scalar field with potential V{{p) = can be described by a perturbed model 

when m <^ 1.) In this case Eq. ( |4.3| ) defines the region of the minisuperspace where the 
perturbative approximation is well-defined. 

Using the canonical transformation 

go = In a, pQ = apa, (4.4) 

the Hamiltonian can be cast in a useful form. We have 

H = ^{pl~pl) + Xe'^^V{q,), (4.5) 

where gi = ip, pi = p^, and V{qi) = V{qi)/\. So the system in the "barred" variables is 
equivalent to a two-dimensional Klein-Gordon particle with a time-dependent potential. In 
these variables the conditions ([4.3|) read 

|A|e6'^o\/(gi) <pg ^p?. (4.6) 

We shall first discuss the system in the case of constant potential. In this case the system 
is integrable and the exact reduction to the Shanmugadhasan form can be found. Then we 
shall move to the case of a potential of the form V^(gi) = g^/n! and apply the perturbative 
techniques. We shall see that for n = the perturbative result coincides with the exact one. 



4.1 Constant Potential: V(qi) = 1 

Let us reduce the system to the Shanmugadhasan form using the method of Sect. |^. The 
generator of the canonical transformation is 



W{q^, P^) = gi^Pi + Po + ^^JPi -Po\Jl + p'^^pf ^° - arctanh ^ 1 + -^^^ee?" 

(4.7) 

Using Eq. ( [4.7| ) we find 



Q = \^± , ^ arcsinh \ ^e-^^o _ 1 , (4.8) 

2 Pi Qpji _ |e65o V2A 
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where different signs refer to "0" and "1" variables respectively. In the next section we 
shall investigate the case of non-constant potential in the limit of small A and compare the 
perturbative result to the constant potential case. So it is convenient to study the behavior 



of Qq in the limit of small A. It is easy to see that the Shanmugadhasan variables ( |4.8| ) 
are ill defined for A ^ 0. This means that this set of variables cannot be used to describe 
the evolution of the system when the potential is close to zero. Since the Shanmugadha- 
san variables are not unique we can redefine them and eliminate the singularity at A = 0. 
Recalling that the generator of the canonical transformation to Shanmugadhasan variables 
is determined up to a generic function of the new momenta we define the new generator 



W\q^,P^)=W{q^,P^)- J dPo 



1 



12 - Po 



In 



Po 



(4.10) 



With this definition the singularity at A = in the Shanmugadhasan variables is removed. 
Pfj, are clearly unaffected. The conjugate variables read now 



1 



In the limit of small A we have 



-3go + In 



1 + . 1 



2A 



2, 1 



(4.11) 



1 ^ go 

2 \Pi Po, 



6 



+ 0{\'). 



(4.12) 



The quantity Qo in Eqs. (|4.11|) can be used as time parameter for the system because 
[Qq,H]p = 1. The Poisson bracket of the first order approximation ( ^.12| ) with the Hamil- 
tonian is 



gl25o 
Po 



(4.13) 



The region of the phase space where the first order approximation for the time parameter 
holds is then defined by the relation 



9A^ 



-4 

Po 



> 



pI 



(4.14) 



or, using the original variables, 

3|A|a^<p^ (4.15) 

The region of validity of the first order approximation is consistent with Eqs. ( [4.3| , [4.6| ) there- 
fore the perturbation techniques can be succesfuUy implemented. In the next subsection we 
discuss the perturbative technique in the case of a nonconstant potential V{qi) = q^/n\ and 
recover the perturbative result of this section in the case n = 0. 
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4.2 Nonconstant Potential: V(qi) = q"/n! 

Let us apply the formalism of Sect. ^ to the case V{qi) = q^/n\. Our starting point is the 
reduction of the Hamiltonian (f4.5|) to the form ( p.2|) . We define the canonical tranformation 
to the new variables 



= iqiTqo)VPiTPo, 
= VPi T Po ■ 

Using Eqs. ( |4.16|J4.1^) the Hamiltonian is cast into the form 

H = po + Ae6(^i-^«)v^^^V [(gi + qo)VlhT]^ 

The first order Hamilton characteristic function can be easily calculated from Eq. 
result is 



W{q„P,) = q,P,- 



X 7" 



+ 0{\' 



(4.16) 
(4.17) 



(4.18) 
. The 

(4.19) 



where f3 = \/Pi — Pq, 7 = \/Pi + Pq, and ^ = qi + qo- Starting from Eq. ( |4.19| ) the first 
order gauge invariant quantities and the first order time parameter can be calculated. We 
have 



Qo = go + 



A 



1 gl2/3gi 



Qi = qi 



2(n- 1)! (-6)"+! 
A 



I gl2/3gi 



2{n - 1)! (-6)"+! 
A 



I2gi7' 


gn 




d(3 




gn 


nl3 ) 


dp 



„2 fin+1 



+ 



(9/5" n/5a/3"+\ 

^2 gn+l 

(9/3" ^ ^9/3"+i 



Po=Po + ^(70"e^'3(«i-««) + 0(A' 
n! 



n! 



n! (-6)" dl3" \ 13 



+ 0(A^ 



+ 0(A^), 
(4.20) 

+ 0(A2), 
(4.21) 
(4.22) 

(4.23) 



where now [3 and 7 are evaluated to the order 0(A), i.e. (3 = y/pi — po, 7 = y/pi + po- In 
the barred variables the above canonical quantities reads 



Qo 
Qi 



A 



\ PI PO J 



1 / gi _ go \ 

2[p, Po) 2(n-l)!(-6)"+i'' 



1 / gi ^ go 



2\pi poj 2(n- 1)! (-6)"+! 



4-^+0(A^) 
4+)+0(A2) 



Po = l{pi'Pl) + ^qie'^'+0{X'), 



Pi = \{p\ + pl) + ^€e' 



2A POP? M'^+'t) 
n\ (-6)" 



5„, + O(A') 



(4.24) 

(4.25) 
(4.26) 
(4.27) 
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where 



Br, 



6pl 



^ + ^ ±1 



/g-6poqi/pi \ 



dpQ 



+ 



p\ d 



,n+l 



npQ dp^^^ 



-6poqi/pi 



Po 



After some algebra and using the Poisson brackets 

Bn,pl ^ = 0, 



Bn^pl 




-12 




P 


V Pi 


' % 
-Dn, — 






Po. 


P 


Po 


-Dn, — 




0, 


Pi 


P 





-6po(ji/pi 



(4.28) 
(4.29) 

(4.30) 
(4.31) 

(4.32) 

(4.33) 



it is easy to verify that [Q/^, Pu]p = ^iiu + O(A^). Hence, Qo defined in Eq. ( [4. 241) is the first 
order approximation of the global time parameter of the system and the quantities defined 
in Eqs. ( [4.25| - ^^27|) form a complete set of approximate gauge invariant quantities. Setting 
n = it is straightforward to check that Eqs. ( |4.24[ - ^^ ) reduce to Eqs. ( [4.12| ) and Eqs. 
(|4.9| ). Therefore the first-order expansions (|4.24| - |4l27|) approximate the exact expression for 
a constant potential. When the potential is not constant the result ( [4.24[py27| ) is of particular 
relevance. In this case the system is not integrable and a global time and gauge invariant 
quantities cannot be found in an analitic form. However, finite expressions like ( [4.24[ - ^^27| ) 
can be found at any order in the A expansion. 

The procedure illustrated above can be easily generalized to the case of a closed FRW 
universe. In this case the Hamiltonian is defined as in Eq. (|4.2| ) with an extra term — a^/2. 
Using the canonical transformation 



1 



In- 



Pa 



pl + a? 



Po = ax/pl + 



(4.34) 
(4.35) 



the kinetic part becomes equal to the kinetic part of Eq. ( |4.5|) and the (perturbed) potential 
part reads 



2po 



e^«"°y(gi) . 



(4.36) 



1 + e^^^o , 

The discussion of this case proceeds along the lines of the flat FRW case, the only difference 
being the form of the perturbed Hamiltonian in the barred variables. 
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5 Conclusions 



In this paper we have illustrated a procedure that allows the approximate quantization of 
finite-dimensional, time-reparametrization invariant, non integrable systems. The method 
is based on the construction of a set of Shanmugadhasan variables approximated to the 
order 0(A") in the perturbation parameter A which identifies the nonintegrable sector of the 
theory. Using the new canonical variables the system becomes trivial and quantization can 
be carried out by standard techniques. 

The quantization scheme illustrated above has some similarities with the method de- 
veloped by Barvinsky [rH| and Barvinsky and Krykhtin |T^. In Ref. [T^ Barvinsky and 



Krykhtin prove the equivalence of BFV, Dirac, and reduced quantization at the one-loop 
(semiclassical) approximation by the equivalence of the physical inner product in the two 
approaches. In the present paper we have proved the equivalence of Dirac and reduced 
methods by the same observation (see below Eq. (p.l8|) ). However, in our case the result is 
obtained by considering not a semiclassical expansion, but a perturbative expansion around 
the integrable sector of the theory. Further, the equivalence of Dirac and reduced approaches 
is proved to an arbitrary order in the perturbation parameter. A deeper discussion on the 
relation between our method and the approach illustrated in Ref. |jl9| is certainly worth 
being investigated. 

What are advantages and disadvantages of the quantization scheme presented in the 
previous sections? We have seen that Dirac and reduced quantization coincide - to the order 
0(A"') - when (approximate) Shanmugadhasan variables are used. Further, the observation 
that Shanmugadhasan variables represent a preferred choice of coordinates in the phase 
space may solve, at least partially, the ambiguity in the choice of gauge fixing and the 
problem of inequivalence of quantum theories obtained starting from different classical sets of 
canonical coordinates. In a sense, Shanmugadhasan variables play the role of the "Cartesian 
coordinates" in elementary quantum mechanics: quantization prescriptions must be given in 
the Shanmugadhasan representation. One of the main limitations of the method is that the 
relation between Shanmugadhasan and original canonical variables is generally nonlinear. 
As a consequence, it is very difficult, if not impossible, to write down quantum operators 
corresponding to the original fields. However, let us stress that one should be very careful 
in assigning a privileged meaning to the original fields. In general relativity no one can 
identify a priori which quantities determine the physical properties of the system: it is the 
nature of the system itself that selects "the physical" fields. In the canonical description 
of the Schwarzschild black hole, for instance, the relevant field is not the metric tensor but 
the mass of the black hole (see Ref. |]2D[), i.e a gauge invariant, nonlinear, function of the 
"original field" (the metric) that appears in the Einstein-Hilbert action. 
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